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Abstract 

Background: Iron oxy hydroxides are commonly found in natural aqueous systems as nanoscale particles, where 
they can act as effective sorbents for dissolved metals due to their natural surface reactivity, small size and high 
surface area. These properties make nanoscale iron oxyhydroxides a relevant option for the remediation of water 
supplies contaminated with dissolved metals. However, natural geochemical processes, such as changes in ionic 
strength, pH, and temperature, can cause these particles to aggregate, thus affecting their sorption capabilities and 
remediation potential. Other environmental parameters such as increasing salinity may also impact metal retention, 
e.g. when particles are transported from freshwater to seawater. 

Results: After using synthetic iron oxyhydroxide nanoparticles and nanoparticle aggregates in batch Zn(ll) adsorption 
experiments, the addition of increasing concentrations of chloride (from 0.1 M to 0.6 M) appears to initially reduce Zn 
(II) retention, likely due to the desorption of outer-sphere zinc surface complexes and subsequent formation of aqueous 
Zn-CI complexes, before then promoting Zn(ll) retention, possibly through the formation of ternary surface complexes 
(supported by EXAFS spectroscopy) which stabilize zinc on the surface of the nanoparticles/aggregates. In batch Cu(ll) 
adsorption experiments, Cu(ll) retention reaches a maximum at 0.4 M chloride. Copper-chloride surface complexes are 
not indicated by EXAFS spectroscopy, but there is an increase in the formation of stable aqueous copper-chloride 
complexes as chloride concentration rises (with CuCI + becoming dominant in solution at -0.5 M chloride) that would 
potentially inhibit further sorption or encourage desorption. Instead, the presence of bidentate edge-sharing and 
monodentate corner-sharing complexes is supported by EXAFS spectroscopy. Increasing chloride concentration has 
more of an impact on zinc retention than the mechanism of nanoparticle aggregation, whereas aggregation condition 
is a stronger determinant of copper retention. 

Conclusions: Based on these model uptake/retention studies, iron oxyhydroxide nanoparticles show potential as a 
strategy to remediate zinc-contaminated waters that migrate towards the ocean. Copper retention, in contrast, appears 
to be optimized at an intermediate salinity consistent with brackish water, and therefore may release considerable 
fractions of retained copper at higher (e.g. seawater) salinity levels. 
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Background 

Anthropogenic activities such as metal ore mining are 
known to cause metal contamination of the surrounding 
environment. Precipitation events and regular seepage/ 
flow in these areas mobilize acidic water and dissolved 
metals, contaminating nearby water sources [1,2]. While 
the harmful effects of some metals such as mercury and 
lead are well known, metals that are common dietary 
requirements such as zinc and copper can also cause 
health problems upon elevated exposure. For example, 
excess copper intake can induce neurological and psy- 
chological problems, while elevated zinc can impair the 
absorption of other ions such as copper and iron and 
cause corrosive damage to soft tissues [3]. Metal con- 
tamination can also disrupt the ecosystems in proximity 
of these water supplies; for example, in plants excess 
zinc is found to inhibit many metabolic functions while el- 
evated copper levels induce injury and oxidative stress, also 
disturbing metabolism and macromolecular activity [4]. 

Iron oxyhydr oxides, specifically goethite (a-FeOOH) 
and ferrihydrite ((Fe) 2 O 3 »0.5 H 2 0), form readily in acid 
mine drainage environments. In these systems, iron oxy- 
hydroxides form when pyrite reacts upon exposure to 
oxygen and water [5], producing large volumes of acid 
that further facilitate the dissolution of sulfides and the 
mobilization of trace metals. Iron oxyhydroxides occur 
at the nanoscale in natural environments [6-8] and are 
inherently effective sorbents for dissolved metals because 
of their small size, high surface area, and natural surface 
reactivity [9-13]. Iron oxyhydroxides are also considered 
to be the predominant reactive mineral phase in lake 
and marine sediments [8], where they can play a signifi- 
cant role in natural attenuation processes. 

Typically, iron oxyhydroxide nanoparticles aggregate 
rapidly [14,15] upon exposure to natural geochemical con- 
ditions based on ionic strength, pH, and temperature, al- 
though the mechanisms of aggregation vary for each of 
these parameters. Increasing pH decreases the particles' 
surface charge density as pH approaches the point of zero 
charge, which for iron oxyhydroxides is between pH 7.0 
and 9.0 [16]. As the surface charge decreases, electrostatic 
repulsion between particles is reduced, allowing them to 
move closer together and aggregate [17-19]. Alternatively, 
increasing ionic strength decreases the electrical double 
layer thickness needed to offset the particles' surface 
charge [20-22], enabling particle aggregation. Lastly, in- 
creasing temperature increases thermal motion of the 
nanoparticles, causing them to collide more often and 
with more force, resulting in aggregation [23] . 

These differences in aggregation mechanism lead to 
the formation of morphologically distinct aggregates. For 
example, increased pH and ionic strength induce rapid 
aggregation, leading to the formation of disordered fractal 
aggregates [23,24], while increased temperature allows 



aggregation to occur more slowly, resulting in the forma- 
tion of more ordered and oriented aggregates [23]. The 
differing morphologies of nanoparticle aggregates exposed 
to varying geochemical conditions can affect their metal 
retention properties by reducing available reactive surface 
area, inhibiting (or creating) access to aggregate interpore 
spaces, and by altering the proportions of different surface 
binding sites [25-27]. 

The salinity of the aqueous environment can also 
affect metal adsorption/retention capabilities. Although 
it has been shown that the presence of chloride typically 
inhibits metal uptake onto iron oxyhydroxides through 
the formation of stable metal-chloride aqueous species 
[12,20,28], less is known about the retention behavior of 
metals that were initially sorbed to nanoparticles or their 
aggregates in a freshwater environment before then be- 
ing exposed to increasing salinity. Simulating the transi- 
tion from freshwater streams and rivers (represented by 
0 M chloride) to brackish (mixed) waters of bays and estu- 
aries (0.1-0.4 M chloride) to marine waters (0.6 M chlor- 
ide) will enable a better understanding of how increasing 
salinity levels affect metal retention to the nanoparticles. 
The effect of nanoparticle aggregation, representing the 
conditions that the particles were exposed to prior to 
metal uptake, is an additional determinant potentially con- 
trolling metal adsorption and retention. 

Previous studies indicate that different aggregation 
mechanisms produce nanoparticle aggregates with vary- 
ing metal retention capabilities [23,29]. This study will 
investigate the effect of chloride concentration and a 
subsequent lowering of pH on Zn(II) and Cu(II) reten- 
tion to unaggregated iron oxyhydroxide nanoparticles as 
well as those aggregated under conditions of elevated 
ionic strength, pH, and temperature. Our hypothesis is 
that introducing metal-sorbed iron oxyhydroxide nano- 
particles or nanoparticle aggregates to an increasingly 
saline environment, as in the transition from freshwater 
to seawater, will systematically reduce Zn(II) and Cu(II) 
retention onto nanoparticle surfaces. The expectation is 
also that more aggregated particles will retain less metal 
initially but will also be less influenced by increasing sal- 
inity, thus exhibiting greater retention than unaggre- 
gated or less aggregated nanoparticles. 

Experimental methods 

Iron oxyhydroxide nanoparticle synthesis 

Iron oxyhydroxide nanoparticles were prepared using a 
flash microwave synthesis technique described by Guyodo 
et al. [30]. Equal volumes of 0.20 M Fe(N0 3 ) 2 and 0.25 M 
NaHC0 3 solutions were prepared before adding the 
NaHC0 3 to the Fe(N0 3 ) 2 dropwise through a 0.20 [im 
syringe filter. Once mixing was complete, the solution 
was agitated on a shaking table and vented periodically 
until the newly formed C0 2 in solution was released 
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(approximately 5 minutes total). After degassing, the so- 
lution was heated in a conventional microwave at high 
intensity for approximately 3.5 minutes in 30-second in- 
tervals just until the onset of boiling to induce nucle- 
ation. To halt nucleation, the resulting nanoparticle 
suspension was immediately placed in an ice bath until 
it had returned to room temperature. The suspension 
was transferred into 1000 MWCO dialysis tubing and 
allowed to equilibrate in deionized water which was re- 
placed 3 times/day until the pH and conductivity had sta- 
bilized at 5.0 and 1.5 (iS/cm, respectively (~5 days). 
Following equilibration, the resulting nanoparticle suspen- 
sion (final solids concentration: 6.7 g/L) was refrigerated 
at 4°C in sealed HDPE bottles until used in aggregation 
and/or batch uptake experiments. 

Nanoparticle aggregation 

Separate aliquots of the initial nanoparticle suspension 
were exposed to varying pH, ionic strength, and tempera- 
ture conditions in order to induce aggregation (Table 1). 
To aggregate the particle suspension under increased pH 
and ionic strength conditions, it was first transferred into 
lengths of 1000 MWCO dialysis tubing and placed into a 
control solution (0.001 M NaNO s and pH 5.0) for 3 days 
to allow equilibration. The particle suspensions, still in the 
same dialysis tubing sections, were then transferred into 
separate solutions of their respective aggregation condi- 
tions (Table 1). The aggregation solutions were replaced 
daily. After 5 days, the particle suspensions were returned 
to control solution to allow re-equilibration for 3 days, 
with the control solution also replaced daily. The particle 
suspensions were then refrigerated at 4°C until they were 
characterized and used in uptake experiments. 

To aggregate the particle suspensions at elevated tem- 
peratures, aliquots of the suspension were placed into 
tightly-capped HDPE bottles. The bottles were placed in 
ovens at their respective temperatures (Table 1) for 
4 days. After 4 days of heating, the particle suspensions 
were refrigerated at 4°C until their use in further 
experiments. 

Table 1 Experimental aggregation conditions listing ionic 
strength, pH, and temperature 



Aggregation condition NaN0 3 concentration pH Temperature 



Control 


0.001 M 


5.0 


RT (~20°C) 


pH 8 


0.001 M 


8.0 


RT (~20°C) 


pH 10 


0.001 M 


10.0 


RT (~20°C) 


0.1 M 


0.1 M 


5.0 


RT (~20°C) 


1.0 M 


1.0 M 


5.0 


RT (~20°C) 


25°C 


0.001 M 


5.0 


25°C 


50°C 


0.001 M 


5.0 


50°C 


75°C 


0.001 M 


5.0 


75°C 



Characterization of nanoparticle aggregates 
X-ray diffraction 

Powder X-ray diffraction patterns were collected on con- 
trol (unaggregated) nanoparticles as well as those aggre- 
gated at 0.1 M ionic strength, 1.0 M ionic strength, 
pH 8.0, and pH 10.0 to assess the crystallinity and min- 
eral phase of the samples. Samples were air-dried in 
50 mL Falcon centrifuge tubes and ground with an agate 
mortar and pestle prior to loading as a thin film between 
two layers of Scotch tape. X-ray diffraction patterns were 
collected at the Stanford Synchrotron Radiation Light- 
source (SSRL) on beamline 11-3 at an energy of 12735 eV 
using a Si(311) crystal monochromator calibrated with a 
powdered LaB 6 standard. Data was collected using a 
Mar345 CCD detector for 90 seconds for all samples ex- 
cept for the pH 8.0 sample, for which data was collected 
for 30 seconds and scaled accordingly. The resulting pat- 
terns were analyzed with the program fit2D [31] and 
background-subtracted using the program XRD-BS [32]. 

Dynamic light scattering 

Dynamic light scattering analysis was conducted on the 
control nanoparticles and all aggregates to examine and 
compare their hydrodynamic diameters. Aliquots of the 
particle suspensions were diluted lOx with deionized 
water for optimal particle detection and placed into 
cuvettes, which were loaded directly into a Malvern 
Zetasizer Nano S dynamic light scattering unit. The 
samples were agitated directly in their cuvettes multiple 
times with a 1000 [iL pipet immediately before analysis 
in an effort to minimize the effect of particle settling. 
Each sample was analyzed between 2-6 times with 70 
cumulative measurements collected per trial, with the 
entire analysis taking approximately 1 minute. 

Geochemical modeling 

Geochemical modeling was performed using thermo- 
dynamic equilibrium constants for Zn(II) [33] and Cu(II) 
[34], producing speciation diagrams representing the pa- 
rameters of the experiment. The concentrations of aque- 
ous Cu(II) and Zn(II) used were based upon the average 
percent uptake following the adsorption phases of the 
macroscopic experiments (0.131 mM Cu(II) and 0.046 mM 
Zn(II) left in solution). The fraction of each metal- 
chloride compound was calculated as a function of 
environmental chloride concentration using the rele- 
vant thermodynamic equilibrium equations and p values, 
finding the concentration of each compound by factor- 
ing the solubility product of sodium chloride into the 
equilibrium. 

Metal adsorption/desorption studies 

Prior to initiating metal uptake, the iron concentrations of 
the unaggregated control suspension and each aggregate 
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nanoparticle suspension were measured as a proxy for 
nanoparticle concentration with a Thermo Scientific 
SOLAAR M Series atomic absorption spectrometer. The 
particle suspension was diluted in a 1:216 ratio in water 
acidified with ultrapure nitric acid (pH <2.0, 0.03 M) in 
order for the sample iron concentrations to fall within the 
optimal (linear) detection range of the AA spectrometer. 
Iron concentrations of the aliquots of aggregate suspen- 
sions ranged from 3365 to 3654 ppm. This mea- 
surement captures sampling variability as a result of 
particle aggregation, allowing for normalization based 
on iron concentration in each aliquot of aggregate sus- 
pension added to the experimental setup. The amount 
of each suspension added to the experimental setups 
was correspondingly adjusted in order to deliver a con- 
sistent quantity of nanoparticles to each reaction vessel, 
minimizing the likelihood that differences in uptake and 
desorption behavior between the various aggregates were 
a function of varying nanoparticle concentration. 

Once normalized, iron oxyhydroxide nanoparticle sus- 
pensions were exposed to a 5 mM metal (Cu(II) or Zn 
(II)) stock solution, adding appropriate volumes of DI 
water, a stock 0.1 M NaN0 3 solution, and 0.1 M NaOH 
and/or 0.1 M HN0 3 to achieve a final target volume of 
150 mL, metal concentration of 0.5 mM, and nanoparti- 
cle solids concentration of 0.17 g/L while also maintain- 
ing ionic strength and pH at control conditions of 
0.001 M and pH 5.0. All experiments were conducted in 
triplicate. After introducing the particle suspensions to 
the dissolved metal solution, the pH was raised with 
20 uL aliquots of 0.1 M NaOH to 6.0 ±0.1 for Cu(II) 
samples and 7.0 ±0.1 for Zn(II) samples to allow for 
maximum metal uptake. These pH values were deter- 
mined by preliminary pH-based uptake experiments and 
calculations as well as those of other investigators [16]. 
The samples were then sealed in HDPE bottles and 
placed on a rotating table for 24 hours. 

After the adsorption step, the samples were split into 
three separate 50 mL aliquots. One aliquot was immedi- 
ately analyzed to assess the initial extent of metal uptake, 
referred to from this point forward as the "adsorbed" 
sample. Appropriate amounts of solid NaCl were added 
to the remaining aliquots to raise the salinity to either 
0.1 M, 0.4 M, or 0.6 M ("salt added" samples). For the 
"salt added" samples, the target salinity level was main- 
tained for 24 hours. For the remaining aliquot, after the 
salinity increase the pH was lowered to 5.0 ±0.1 using 
20 uL aliquots of 0.1 M HN0 3 ("pH dropped" samples) 
and the suspension agitated for 24 hours. While kinetics 
were not monitored in this study, the 24-hour exposure 
periods were determined based on our previous ion se- 
lective electrode lab studies which showed that adsorp- 
tion and desorption are each complete after one hour. 
Exposure periods of 24 hours were used to better 



simulate environmental conditions and to ensure thor- 
ough adsorption and desorption. 

Following exposure, samples were centrifuged at 3000 rpm 
in 50 mL Falcon tubes for 10 minutes. The resulting 
supernatant was decanted, filtered using 0.20 um syringe 
filters, diluted lOx in acidified (pH < 2.0) DI water, and an- 
alyzed for either Zn(II) or Cu(II) using AA spectrometry. 
The metal concentration obtained from AA analysis was 
used to calculate the degree of metal uptake/retention to 
the substrates following each of the experimental steps, 
assuming minimal uptake of dissolved metal to the vessel 
walls. 

Extended X-ray absorption fine structure (EXAFS) 
spectroscopic analysis 

Selected solid samples from the Zn(II) and Cu(II) 
macroscopic experiments were collected for EXAFS ana- 
lysis. The solids were spread onto Whatman filter paper 
to remove excess water and loaded into Teflon sample 
holders with Kapton tape as moist pastes for analysis. 

All Zn(II) K-edge EXAFS data were collected on the 
moist pastes at SSRL on beamline 4-3 in fluorescence 
mode at room temperature using a 13-element high- 
throughput germanium detector. Cu(II) K-edge EXAFS 
data were collected at SSRL on beamline 11-2 using a 
100-element high-throughput germanium detector in 
fluorescence mode. The fluorescence method is advanta- 
geous for lower concentration samples [35,36]. Aluminum 
filters were used to reduce Fe K-edge fluorescence and 
zinc and copper metal foils were used as calibrants. 

The resulting spectra were analyzed using the SIXPack 
data processing software version 1.01 [37]. Deadtime 
corrections were performed on each scan in order to ac- 
commodate for loss of signal upon saturation of the de- 
tector channels prior to being averaged together. Those 
averages were converted to k-space with a k 3 weighting 
and were Fourier transformed. All copper average files 
(excluding the spectra from the pH 10 aggregates at 0.4 
M chloride) required minor deglitching due to mono- 
chromator imperfections. Both Zn and Cu EXAFS spec- 
tra were fit over a k- range of 3.0-12.0 A using model 
scattering paths which were generated in SIXPack using 
Fefiftl [38]. 

Results and discussion 

Nanoparticle characterization 
X-ray diffraction 

Generally, the XRD patterns obtained from the aggre- 
gated particles are similar to those of the unaggregated 
control particles (Figure 1), suggesting that exposing the 
nanoparticle suspension to aggregation conditions does 
not significantly alter the mineral phase or degree of crys- 
tallinity. The scattering pattern of the pH 10-aggregated 
particles shows higher peak intensities consistent with 
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10 20 30 40 50 



29 (deg, 12735 eV) 

Figure 1 Powder X-ray diffraction patterns of the control, ionic 
strength and pH aggregated particles. Ferrihydrite (solid gray 
lines [18]) and goethite (dashed lines, No. 29-0713) PDFs are included 
for comparison. 



larger aggregates. X-ray diffraction analysis was not per- 
formed on temperature-aggregated particles in this study; 
however, based on our previous studies [29], aging these 
same iron oxyhydroxide nanoparticles at 78°C caused a 
gradual transformation from ferrihydrite to goethite over 
a 7-day duration. 

Comparison of the collected diffraction patterns to PDF 
standards of 6-line ferrihydrite and goethite (Figure 1) in- 
dicates that the nanoparticle suspensions are primarily 
comprised of 6-line ferrihydrite. However, there are a few 
specific peaks present that are unique to goethite. The 
particles therefore appear to be nanoparticulate ferrihy- 
drite, but some proportion of the sample is goethite or has 
some degree of goethite-like structure. 

Dynamic light scattering 

Analysis of the nanoparticle aggregates with dynamic 
light scattering demonstrates that as the nanoparticles 
are exposed to elevated pH, ionic strength, or tempera- 
ture, their z-diameter increases, indicating the formation 
of larger aggregates (Figure 2). This confirms that aggre- 
gation is induced by these parameters and further veri- 
fies that returning the aggregates to control solution 
following exposure to aggregation conditions is not suffi- 
cient to provoke significant reversal of the aggregation. 

The results also suggest that elevated pH or ionic 
strength conditions produce larger aggregates than ele- 
vated temperature over the timeframes investigated and 
is indicative of different aggregation mechanisms and 




Control 0.1 M1 .0 M pH 8pH 10 25°C 50°C 75°C 



Experimental Condition 

Figure 2 Hydrodynamic diameters of control and aggregated 

nanoparticles measured using dynamic light scattering analysis. 

Results are shown on a log scale, 
k J 

rates. When the nanoparticles are placed in environ- 
ments with elevated temperature, they aggregate slowly 
(and often in an oriented manner), forming relatively 
compact, ordered aggregates. In contrast, nanoparticles 
exposed to elevated pH or ionic strength aggregate more 
rapidly, forming disordered aggregates with substantial 
interstitial water and resulting in a larger aggregate 
diameter [23]. Dynamic light scattering results are con- 
sistent with X-ray diffraction results, both of which sug- 
gest that higher pH, ionic strength, and temperature 
induce greater degrees of aggregation, with particles ag- 
gregated at pH 10 being the largest. Earlier small and 
wide-angle X-ray scattering studies [23] also support 
these findings of increased aggregation with elevated pH, 
ionic strength, and temperature. 

Zinc sorption 
Macroscopic results 

Raising the pH to 7.0 ±0.1 in the adsorption phase 
caused 85-95% of the zinc in solution to sorb to the par- 
ticle surfaces ("adsorbed" samples, Figure 3a-d). On 
average, the initial addition of chloride caused minor de- 
sorption to occur (0.1 M "salt added" samples) relative 
to the substantial decrease in zinc retained when the pH 
was lowered in the presence of chloride ("pH dropped" 
samples). In both the "salt added" and "pH dropped" 
samples, a general increase is seen in the percent zinc 
retained as chloride concentration increases, following a 
slight drop in percent zinc retained upon the addition of 
0.1 M NaCl, with the trend appearing much more con- 
sistently in the "pH dropped samples" (Figure 3a-d). 
These results suggest that increasing chloride concentra- 
tion stabilizes the zinc that is sorbed to the particle sur- 
faces, causing more zinc to be retained. 

The initial decline in zinc retention after lowering pH 
between the salt-free and 0.1 M salinity trials (Figure 3a-d) 
suggests that it may be more thermodynamically favorable 
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Figure 3 Macroscopic zinc uptake/retention data 
d) temperature-aggregated nanoparticles. Adsorbed 
columns) represent pH 6.0 ± 0.1 and 0.1-0.6 M NaCI, and 
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for free or weakly-bound zinc, such as zinc in outer- 
sphere sorption complexes, to form stable aqueous zinc- 
chloride species over surface sorption complexes within 
this salinity range. As the chloride concentration rises, 
however, more chloride is available to react with zinc 
sorbed through inner-sphere mechanisms. It is possible 
that chloride modifies the type of surface complex formed 
between the zinc and the nanoparticle, stabilizing the Zn- 
surface bond. 

Comparing only the retention data from the "pH 
dropped" samples highlights the initial decline in reten- 
tion from the salt- free to 0.1 M Cl~ conditions and the 
progressive increase in Zn(II) retention as the salinity is 
further increased to 0.6 M (Figure 4). Upon introducing 
a saline environment, the differences in retention be- 
tween aggregates decreased (as assessed by the standard 
errors within each salinity category shown in Table 2, 
calculated using results from all samples), indicating 
that the effect of increased environmental chloride 



No salt 



0.1 M CI 



0.4 M CI 



0.6 M CI 
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Figure 4 Macroscopic zinc retention data for all samples at all 
salinities following exposure to NaCI for 24 hours at pH 5.0 ± 0.1 . 

Salinity increases from left to right as labeled. 
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Table 2 Average percent Zn(ll) uptake and standard 
errors for "pH dropped" samples 





Average (% uptake) 


Standard error 


No salt 


39.1% 


3.14% 


0.1 M CI" 


28.1% 


0.96% 


0.4 M CI" 


36.7% 


0.47% 


0.6 M cr 


63.2% 


0.73% 



concentration is stronger than the effect of the aggrega- 
tion method on the retained fraction. 

Geochemical modeling 

As the environmental chloride concentration increases, 
the fraction of Zn 2+ in solution decreases as aqueous zinc- 
chloride complexes are formed (Figure 5). ZnCl + forms 
the fastest, overtaking Zn 2+ as the dominant aqueous zinc 
species at -0.65 M chloride, although ZnCl 2 and ZnCl 3 " 
are also present in smaller proportions (ZnCl 4 2 " is not 
shown because its corresponding p value was not pro- 
vided). This experiment only tested zinc retention up to 
0.6 M chloride, but previous macroscopic studies con- 
ducted in our lab suggest that zinc retention plateaus at 
chloride concentrations between 0.6 M and 1.2 M, which 
could correlate with the dominant zinc species becoming 
ZnCT at 0.65 M chloride. 

Spectroscopic results 

Subtle spectral differences can be observed between the 
adsorbed samples and those whose pH levels were low- 
ered in the presence of chloride (Figure 6). These differ- 
ences are most apparent in the Fourier transforms, in 
which the second-neighbor features appear to merge to- 
gether as chloride concentration increases. The EXAFS 
fitting results of the first-neighbor feature show an 
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Figure 5 Speciation diagram of Zn(ll) chloride species at a 
range of chloride concentrations and a Zn(ll) concentration of 
0.046 mM at pH 5.0. 
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Figure 6 Zinc K-edge EXAFS spectra and Fourier transforms 
(black) with overlain fits (grey) for "adsorbed" and "pH dropped" 
states (0.1 M CI" and 0.6 M CI") of samples aggregated at 1.0 M, 
75°C, and pH 10. 



average coordination number of 3.6 (range: 3.6-3.7) 
(Table 3) for the Zn-O shell in the adsorbed samples, in- 
dicating that in the adsorption phase, sorbed zinc is ini- 
tially binding to the nanoparticle surfaces in a 
dominantly tetrahedral coordination environment. These 
results, in conjunction with previous studies and fitted 
Zn-O bond distances of 1.99 ± 0.01 A, suggest that zinc 
is forming bidentate corner-sharing complexes with the 
iron octahedra that comprise the nanoparticle surfaces 
[29,39], consistent with zinc binding to ferrihydrite 
[39-42]. The increase in the Zn-O coordination number 
of the "pH dropped" samples (average: 5.4, range: 5.2- 
6.1) compared to that of the "adsorbed" samples (aver- 
age: 3.6) also indicates a shift from the zinc being 
bound on less ordered, tetrahedrally-coordinated sites to 
more ordered, octahedrally-coordinated sites. 

The adsorbed samples were best fit with two Zn-Fe 
shells, the first with an average coordination number of 
0.9 (range: 0.9-1.0) and bond length of 3.17 A (range: 
3.16-3.18 A), and the second with an average coordin- 
ation number of 1.1 (range: 1.0-1.1) and bond length of 
3.42 A (range: 3.41-3.43 A) (Table 3). These results gen- 
erally correlate with Juillot et al. [39], who reported bond 
lengths of 3.47 A for zinc bonded to ferrihydrite and 
3.07 A and 3.26 A for goethite. The long bonds present 
in the ferrihydrite samples indicate the presence of 
bidentate corner-sharing complexes [39,43], while the 
shorter Zn-Fe bonds in the goethite samples correspond 
to bidentate (3.26 A) and tridentate (3.07 A) face- 
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Table 3 Results from K-edge EXAFS fitting of Zn(ll) samples (see Figure 6 for corresponding EXAFS spectra and 
Fourier transforms) 



Aggregation condition 


Sorption step 


Zn-0 


Zn-CI 


Zn-Fe 1 


Zn-Fe 2 


R-Factor 


CN 


R (A) 


CN 


R (A) 


CN 


R (A) 


CN 


R (A) 


1.0 M 


Adsorbed 


3.6 ± 0.5 


1.98 ±0.01 






0.9 ± 0.4 


3.1 8 ±0.05 


0.5 ± 0.3 


3.39 ± 0.05 


0.0027 


1.0 M 


0.1 M CI" 


5.2 ±1.0 


2.01 ±0.01 


0.8 ± 0.4 


2.99 ± 0.06 






1 .4 ± 0.5 


3.46 ± 0.04 


0.0045 


1.0 M 


0.6 M cr 


5.4 ±0.8 


2.02 ±0.01 


1 .8 ± 0.5 


2.99 ± 0.05 


1 .6 ± 0.5 


3.05 ± 0.05 






0.0024 


75°C 


Adsorbed 


3.6 ± 0.4 


1.98 ±0.01 






0.9 ± 0.3 


3.1 6 ±0.03 


0.4 ± 0.2 


3.42 ± 0.04 


0.0014 


75°C 


0.1 M cr 


5.2 ±0.9 


2.03 ±0.01 


0.7 ± 0.3 


2.95 ± 0.05 






0.5 ± 0.4 


3.45 ± 0.07 


0.0039 


75°C 


0.6 M cr 


5.2 ±0.8 


2.04 ±0.01 


1 .3 ± 0.5 


3.01 ±0.06 


0.9 ± 0.5 


3.07 ± 0.07 






0.0032 


pH 10 


Adsorbed 


3.7 ±0.3 


1.98 ±0.01 






1 .0 ± 0.3 


3.1 8 ±0.03 


0.6 ± 0.3 


3.39 ± 0.05 


0.0012 


pH 10 


0.1 M cr 


6.1 ±1.2 


2.02 ± 0.02 


1 .0 ± 0.4 


2.97 ± 0.05 






1.1 ±0.5 


3.45 ± 0.04 


0.0041 


pH 10 


0.6 M cr 


5.5 ±0.6 


2.03 ±0.01 


1 .8 ± 0.4 


2.99 ± 0.04 


1 .7 ± 0.4 


3.05 ± 0.03 






0.0015 



Included are coordination numbers (CN), interatomic distances (R), and goodness of fit values (R-factor). Debye-Waller values were allowed to float for first-shell fits 
(average: 0.004 A 2 ) and were fixed at 0.01 A 2 for subsequent shells. 
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Figure 7 Macroscopic copper uptake/retention data for a) control, b) pH-aggregated, c) ionic strength-aggregated, and 

d) temperature-aggregated nanoparticles. "Adsorbed" samples (green columns) represent pH 7.0 ± 0.1, 0.1 M NaN0 3 , "salt added" samples (yellow 

columns) represent pH 7.0 ±0.1 and 0.1-0.6 M NaCI, and "pH dropped" samples (red columns) represent pH 5.0 ±0.1, 0.1 M NaN0 3 , 0.1-0.6 M NaCI. 
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Figure 8 Macroscopic % copper retention data for all particles 
at all salinities examined following exposure to NaCI for 24 hours 
at pH 5.0 ±0.1. Salinity increases from left to right as labeled. 



sharing complexes [39,44]. The presence of both ferrihy- 
drite and goethite binding sites is consistent with the 
X-ray diffraction results, which indicate that the nano- 
particle suspensions exhibit features consistent with both 
ferrihydrite and goethite (Figure 1). The adsorbed sam- 
ples were best fit with bond lengths correlating with past 
studies on both goethite and ferrihydrite, and indicate 
contributions from bidentate edge-sharing (3.17 A) and 
bidentate corner-sharing (3.42 A) complexes (Table 3). 

As the environmental chloride concentration is in- 
creased, there are consistent changes to the Zn-Fe fits. 
At 0.1 M chloride, the best fits resulted in an average 
Zn-Fe coordination number of 1.0 (range: 0.5-1.4) and 
bond length of 3.46 A (range: 3.45-3.46 A) (Table 3), 
which most accurately indicates bidentate corner- sharing 
complexes, as seen in previous ferrihydrite studies [39,43] . 
With increasing chloride concentrations, the observed 
Zn-Fe bond lengths shorten (average: 3.06 A, range: 
3.05-3.07 A) and coordination numbers increase (average: 
1.4, range: 0.9-1.7) (Table 3), corresponding to the presence 
of tridentate face-sharing complexes [39,44]. 

The identification of Zn-Cl neighbors in the EXAFS 
fitting results indicates that zinc is on average bound to 
more chloride ions as the chloride concentration of the 
experimental system increases, with the coordination 
number approaching a value of 2. The Zn-Cl coordin- 
ation numbers indicate that octahedrally-coordinated 
zinc has up to 2 chloride ligands, which would reduce 

Table 4 Average percent Cu(ll) uptake and standard 
errors for "pH dropped" samples 

Average (% uptake) Standard error 



No salt 
0.1 M CI" 
0.4 M CI" 
0.6 M CI" 



60.4% 
57.8% 
71.1% 
55.7% 



1 .54% 
2.19% 
2.34% 
3.88% 
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Figure 9 Speciation diagram of Cu(ll) chloride species at a 
range of chloride concentrations and a Cu(ll) concentration of 
0.131 mM and pH 5.0. 



surface charge repulsion and sorb as an uncharged tern- 
ary surface complex. An average Zn-Cl interatomic dis- 
tance of 2.98 A (2.95-2.99 A) is consistent with a direct 
Zn-Cl bond [45], supporting the conclusion that zinc 
and chloride are binding directly to each other (Table 3). 
The results from the macroscopic experiments and 
EXAFS analysis suggest that the octahedral surface com- 
plexes bound to more ordered binding sites showing evi- 
dence of a direct Zn-Cl bond and tridentate face-sharing 




3579 11 01234 

k(A 1 ) R + A (A) 

Figure 10 Copper K-edge EXAFS spectra and Fourier transforms 
(black) with overlain fits (grey) for "adsorbed" and "pH dropped" 
samples (0.4 M CI" ) for particles aggregated at 1.0 M, 75°C, 
and pH 10. 
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Table 5 Results from K-edge EXAFS fitting of Cu samples (see Figure 10 for corresponding EXAFS spectra and 
Fourier transforms) 

Cu-0 Cu-Fe 1 Cu-Fe 2 

Aggregation condition Sorption step CN R(A) CN R(A) CN R(A) r-factor 



1.0 M 


Adsorbed 


4.3 ± 0.9 


1 .96 ± 0.02 


1.0 M 


0.4 M CI" 


4.7 ± 1 .0 


1 .97 ± 0.02 


75°C 


Adsorbed 


4.4 ± 1 .0 


1 .95 ± 0.02 


75°C 


0.4 M CI" 


4.7 ± 1 .0 


1 .97 ± 0.02 


pH 10 


Adsorbed 


4.5 ± 0.9 


1 .96 ± 0.02 


pH 10 


0.4 M cr 


3.6 ± 0.6 


1.97 ±0.01 



Included are coordination numbers (CN), interatomic distances (R), and goodness 
(average: 0.004 A 2 ) and were fixed at 0.01 A 2 for subsequent shells. 



1 .3 ± 0.6 


2.96 ± 0.04 


0.8 


±0.6 


3.44 


±0.07 


0.0058 


1 .4 ± 0.6 


2.97 ± 0.04 


1.3 


±0.6 


3.44 


±0.05 


0.0059 


1 .4 ± 0.6 


2.94 ± 0.04 


1.3 


±0.7 


3.43 


±0.05 


0.0059 


1 .9 ± 0.6 


2.99 ± 0.03 


1.3 


±0.6 


3.48 


±0.05 


0.0052 


1 .3 ± 0.5 


2.95 ± 0.03 


0.9 


±0.4 


3.44 


±0.06 


0.0052 


1 .5 ± 0.3 


2.99 ± 0.02 


0.9 


±0.3 


3.46 


±0.04 


0.0035 



of fit values (R-factor). Debye-Waller values were allowed to float for first-shell fits 



complexes are more stable, because the conditions under 
which these complexes are formed showed higher per- 
centages of retained zinc. 

Copper sorption 
Macroscopic results 

Following the initial adsorption phase, 70-80% of the 
copper in solution was retained (Figure 7a-d). The sub- 
sequent addition of sodium chloride improves copper re- 
tention for all sets of aggregates. As the salinity was 
increased, the percent copper retained generally reached 
a maximum at 0.4 M chloride. Upon comparing only the 
"pH dropped" samples, the trend showing a retention 
maximum at 0.4 M CI" is more readily apparent 
(Figure 8). In contrast with the Zn(II) results, the differ- 
ences in retention between the various aggregates in- 
crease as salinity rises, as shown by an increasing 
standard error between samples as the chloride concen- 
tration increases (Table 4). These results indicate a fun- 
damentally different behavior of copper compared to the 
zinc experimental results: as chloride concentration in- 
creases, characteristics of the individual nanoparticle 
aggregates, such as morphological, structural, surface 
charge, and surface area differences have a noticeable ef- 
fect on copper retention. 

Geochemical modeling 

While Cu 2+ is the dominant aqueous copper species at 
low chloride concentrations, CuQ + becomes the dominant 
aqueous complex at -0.5 M chloride. In the macroscopic 
results, a retention maximum at 0.4 M chloride was evi- 
dent before a decline in copper retention at 0.6 M chlor- 
ide, which correlates with CuCl + becoming dominant 
in solution at 0.5 M chloride and the significant for- 
mation of other copper-chloride complexes (CuCl 2 , 
CuCl 3 ", and CuCl 4 2 ") (Figure 9). The presence of these 
stable complexes likely hinders further copper sorption, 
leading to the decline in copper retention beyond 0.4 M 
chloride. 



Spectroscopic results 

The Cu K-edge EXAFS spectra of the adsorbed and des- 
orbed samples display visible differences (Figure 10) in- 
cluding a shoulder feature at k = -7.5 A" 1 that is more 
pronounced in the spectra of the desorbed samples. The 
first-shell neighbor was best fit with Cu-O scattering in- 
teractions, with an average coordination number of 4.3 
(range: 3.6-4.7) and an average interatomic distance of 
1.96 A (range: 1.95-1.97 A) (Table 5) for all samples, 
consistent with the equatorial Cu-O bonds in a Jahn-Teller 
distorted octahedral Cu(II) complex [46,47]. These results 
correspond with previous spectroscopy studies of Cu(II) 
sorption onto mineral surfaces [46,48]. 

During EXAFS fitting, the inclusion of a second shell 
Cu-Cl bond with a coordination number of 1 and bond 
length of approximately 2.2 A [49] was attempted to ex- 
plore the possibility of Cu-Cl sorption species. However, 
this caused Cu-O coordination numbers to decline to 
unrealistic numbers and provided unreasonable fits, so 
we have not included Cu-Cl interactions between the 
first Cu-O shell and the second Cu-Fe shell. The second 
shell was best fit with a Cu-Fe neighbor with an average 
coordination number of 1.5 (range: 1.3-1.9) and an aver- 
age interatomic length of 2.97 A (range: 2.94-2.99 A); 
the third shell was also best fit with a Cu-Fe neighbor 
with an average coordination number of 1.3 (range: 0.8- 
1.3) and an average interatomic length of 3.44 A (range: 
3.43-3.48 A) (Table 5). The shorter Cu-Fe distance 
(2.94-2.99 A) has been associated with the formation of 
inner-sphere edge sharing complexes [50]. Similar bond 
lengths have been reported for Cu-Cu bonds corre- 
sponding to the formation of dimers [51], but our best 
fit was obtained with the inclusion of a Cu-Fe shell at 
this length. The longer Cu-Fe distance (3.43-3.48 A) 
most accurately corresponds with the presence of mono- 
nuclear monodentate complexes, as modeled or inter- 
preted by other investigators [52-54]. Visual changes in 
the EXAFS spectra therefore likely correspond to changes 
in the proportions of these species as a result of the 
pH-lowering desorption step. Based on the lack of 
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chloride in the EXAFS fits, it is not likely that chloride 
plays a direct role in the sorption mechanism, e.g. through 
the formation of ternary surface complexes; however, it 
may play an indirect role by initially reducing positive 
surface charges (thereby improving copper retention) 
and then, at higher concentrations, by forming stable 
aqueous Cu-Cl complexes (reducing copper retention). 
Our copper speciation diagram indicates the formation 
of aqueous copper chloride species as chloride concentra- 
tion increases, supporting the relationship between aqueous 
copper speciation and retention behavior (Figure 9). 

Conclusions 

Based on these model studies, iron oxyhydroxide nanopar- 
tides could be a useful tool for removing zinc from con- 
taminated water supplies that lead to the ocean. Copper 
retention, in contrast, appears to be optimized at an inter- 
mediate salinity consistent with brackish water, and there- 
fore may release considerable fractions of retained copper 
at higher (e.g. seawater) salinity levels; copper retention 
also appears to become more variable and dependent on 
aggregation mechanism at these increasing salinities. Ac- 
cordingly, strategies for the environmental remediation of 
metal-contaminated waters should take into account po- 
tential changes in geochemical parameters that may in- 
duce aggregation, increase salinity, and affect solution pH. 
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